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C. T. H. Davies a 

a Department of Physics and Astronomy, University of Glasgow, 
Glasgow G12 8QQ, UK 

Calculations of the heavy-heavy spectrum present a good opportunity for precision tests of QCD using lattice 
techniques. All methods make use of a non-relativistic expansion of the action and its systematic improvement 
to remove lattice artefacts. There was convincing demonstration this year that these methods work and that the 
associated perturbation theory is well-behaved. Comparison to experimental results yields an accurate value for 
the lattice spacing, a, a key result in the determination of a s , and (for the first time this year) the mass of the b 
quark (4.7(1) GeV). 



1. INTRODUCTION 

The spectrum of heavyonium states is an ideal 
place to provide a precision test of QCD from first 
principles. Very accurate calculations are possi- 
ble using lattice techniques. The reason is that 
both statistical and systematic errors can be made 
small, something we are nowhere near achieving 
for light hadron calculations. 

The key ingredients are 

• A non-relativistic formulation appropriate 
to the physics. 

• Moderately-sized lattices at moderate val- 
ues of [3 with perturbatively improved ac- 
tions. 

• High statistics from 0(100) configurations 
and the use of multi-exponential fits to mul- 
tiple correlation functions. 

Results this year demonstrate clearly that per- 
turbatively improved actions work and give con- 
tinuum results when values of the lattice spacing 
are not very small. This is particularly useful for 
heavy quark systems where lattice discretisation 
errors are relatively large. For light hadron sys- 
tems, however, its use would allow calculations 
at lower values of (3 and smaller volumes than are 
currently employed, which might be an advan- 
tage. Another clear result is the greatly increased 
confidence in the ground state that is obtained 
from a multiple exponential fit to multiple corre- 
lation functions. Such a fit is not often possible 



at the level of statistics employed in current light 
hadron calculations. 

Results were presented by four groups this year. 
The groups differ in the action employed for the 
heavy quarks, the level of improvement of the ac- 
tion used, and the gauge field configurations used. 
I shall refer to the collaborations as FNAL 
KEK §, NRQCD ||L UKQCD §,§. 

The FNAL group [|l],|| use improved heavy 
Wilson fermions for both cc and bb spectra on 
quenched gauge field configurations at several val- 
ues of [3, but principally using 100 configurations 
at (3 =6.1. 

The KEK || group use Wilson fermions for 
the cc spectrum but, for the first time, on gauge 
field configurations with 2 flavours of dynamical 
quarks. They have 75 configurations at f3 = 5.7. 

The NRQCD group uses Non-relativistic QCD 
(NRQCD) to study bb § (with 100 quenched 
gauge field configurations at (3 = 6.0 from the 
Staggered Collaboration) and cc spectra Q] (with 
100 quenched gauge field configurations at (3 
= 5.7 from the UKQCD collaboration). The 
NRQCD action is improved at next-to-leading or- 
der. 

The UKQCD collaboration has results for both 
improved heavy Wilson |Q (cc and extrapolating 
to bb) and leading order NRQCD (for bb) || prin- 
cipally at (3 = 6.2 on 60 quenched configurations. 

The most complete spectrum is presented 
by the NRQCD collaboration, with both spin- 
averaged splittings and s and p hyperfine split- 
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tings for bb and cc. I will quote results from other 
groups where available. The KEK group give only 
the 1S-1P spin-averaged splitting which is useful 
for the determination of a -1 and a s . Their re- 
sult provides a confirmation of the unquenching 
corrections that other groups employ to convert 
their quenched a s to a physical one. This is dis- 
cussed by El-Khadra || and I will not refer to the 
KEK results further. 

1.1. NRQCD 

NRQCD makes use of an expansion of the 
quark action in powers of v 2 /c 2 where v/c is a typ- 
ical velocity inside the heavy meson, v 2 /c 2 ~ 0.1 
for T and 0.3 for J/^>. The terms in the action 
can be ordered by power-counting rules || . 

The lowest order terms are simply 



tings. These are given by 



(1) 



if) is a 2-component quark field and the separate 
antiquark has the same action. The presence of a 
single time derivative means that the propagator 
can be calculated in one sweep through the lat- 
tice using a suitable time evolution equation ^ 
derived from a lattice version of the action. 

This leading order action can already give use- 
ful results for spin-averaged splittings, i.e. split- 
tings in which states of a given L are averaged 
over, weighted by the value of 2 J+ 1. These split- 
tings can be compared to those extracted from a 
free Schrodinger equation with the usual central 
heavy quark potential. The leading order action 
has errors of 10% for bb and 30% for cc. The 
UKQCD collaboration work with this action j|] . 

The lowest order spin-dependent terms (those 
with an explicit a are 
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The first of these gives s state hyperfine splittings 
( 3 Si — 1 S'o) and the second p state hyperfine split- 
tings (between x states). Inclusion of these terms 
gives the hyperfine splittings to 10% for bb and 
30% for cc. 

The NRQCD group includes also the next-to- 
leading contributions to the spin-averaged split- 
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(D.E-E.D)ifj (3) 



In addition there are 0(a 2 ) corrections that ap- 
pear at the same level and are included through 
terms containing D A and Df ||. Now spin- 
averaged splittings have systematic errors at the 
1% level for bb and at the 10% level for cc from 
terms neglected in the heavy quark action. 

All the terms above have been written with 
their tree level coefficients. Potentially large ra- 
diative corrections to these coefficients can be ab- 
sorbed by a transformation of the gauge fields : 
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This effectively includes in each term tadpole cor- 
rections to all orders || . These are assumed to be 
the dominant corrections || so that any further 
radiative corrections can be ignored at the next- 
to-leading order level. Explicit calculations |1(J 
for the coefficients of D 4 and Df confirm this. 

1.2. Heavy Wilson 

The approach of heavy Wilson fermions is very 
similar to that of NRQCD (in fact, eventually 
identical). It begins by adding a clover term to 
the usual Wilson action to correct for errors at 
0(a). 



^clover — ~'lrF'ij l i/0'[j,i/'tp 



(5) 



Propagators are calculated by the usual meth- 
ods for Wilson fermions. Convergence for heavy 
quarks (low k) is fast but it still requires O(10) 
passes through the lattice. The FNAL group give 
c the value calculated from the dominant tadpole 
effects described above Qj. This value is Uq = 1.4 
at P = 6.1. The UKQCD group @ give c its tree 
level value of 1.0. 

Compared to the nlo NRQCD approach above, 
the heavy Wilson action has terms missing at 
0(a 2 ) and at order D A . However, there is some 
remnant of the required D A term present so the 
error should be less than the 10 % or 30% of lo 
NRQCD. 
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1.3. Gauge Field Configurations 

There are two sources of systematic error 
present in standard quenched gauge field config- 
urations which must be taken into account when 
assessing the results of the FNAL, NRQCD and 
UKQCD groups. Luckily the two effects tend to 
counteract each other to some extent. 

The first is the presence of 0(a 2 ) errors in the 
usual Wilson plaquette gauge action. This can be 
corrected by the use of 2 x 1 of 2 x 2 plaquettes 
with an appropriate tadploe-improved coefficient. 

We can estimate the size of these effects from 
a Schrodinger equation on the lattice using a lat- 
tice heavy quark potential (such as that from Bali 
and Schilling pJ|). The 0(a 2 ) errors in the gluon 
propagator show up as a lack of rotational invari- 
ance in the potential at small R. If we correct the 
[4sin 2 (ka/2)]~ 1 of the one-gluon exchange term 
to [4sin 2 (ka/2) + 4/3sm 4 (fca/2)] _1 and compare 
the spectra at fixed quark mass, we can obtain 
an estimate of the size of these 0(a 2 ) errors. We 
find a 2% reduction of the spin-averaged IP — IS" 
splitting when the corrections are included. 

Another, more serious, source of error is that of 
the quenching of the gauge fields (at least partly 
removed in the KEK results). We are study- 
ing states mainly well below threshold for de- 
cay to heavy-light channels, so the principal ef- 
fect of quenching is from the change in the shape 
of the heavy quark potential at small distances. 
The quenched coupling constant runs too fast to 
zero, so the potential close to R — is not deep 
enough. Quenching then has the effect of raising 
states which are concentrated at short distances, 
such as s states. A comparison of results from a 
lattice Schrodinger equation with quenched and 
unquenched potentials from the MTc collabora- 
tion |L4| gave a shift to the ratio of 2S — IS* to 
IP — IS" of 10% (downwards on unquenching) . 
The potentials did not have the same lattice spac- 
ing so an absolute shift for \p — IS could not 
be derived. We fixed the quark mass in physical 
units to 4.7 GeV in both cases. 

In perturbation theory the s state hyperfinc 
splitting is proportional to the square of the wave- 
function at the origin and a s (M). In the above 
analysis we did not find a big effect from quench- 
ing on the wavefunction at the origin. This dis- 



agrees with earlier results from El-Khadra (T2J 
who used a continuum Richardson potential. We 
believe that the lattice Schrodinger approach may 
more accurately reflect what we will find in lat- 
tioce simulations. The effect of quenching on 
a s (M) is to reduce it by 20%. Thus the lattice 
Schrodinger approach predicts an s state hyper- 
fine splitting 20% too low in the quenched ap- 
proximation and the continuum Richardson gives 
a 40% reduction. 

2. RESULTS 

There are two parameters to be fixed in all 
these calculations. 

• a -1 : Fix from the spin-averaged spectrum 
which is relatively insensitive to M. (Com- 
pare the s and p splittings for bb and cc from 
the Particle Data Book Q). 

• M : Fix from the kinetic mass, M2, of the 
T or J/^/rj c as appropriate. M2 is given by 
the non-relativistic dispersion relation for 
the energy of a meson at finite momentum: 

m-M 1 + J^- + ... ( 6 ) 

2.1. bb results 

Figure 1 shows the 'spin-averaged' spectrum for 
bb obtained by the NRQCD collaboration at (3 = 
6.0 and the UKQCD collaboration (using leading 
order NRQCD) at (3 — 6.2. Since no experimen- 
tal mass is available for the rjb (IS, or its radial 
excitations, 2S and 3S), no spin-averaging of s 
states is done by NRQCD. The experimental re- 
sults, given as bars, are for the T( 3 Si), IS, 2S 
and 3S. The UKQCD results are inevitably spin- 
averages since they have no spin-dependent terms 
in the action. For the p states a spin-average over 
the 3 Po,i,2 states, all known experimentally at IP 
and 2P, is done. In addition the NRQCD collabo- 
ration have a mass for the as yet unseen 1 P\ state, 
the hb- 

By performing a bootstrap fit to the whole 
spectrum of Figure 1, the NRQCD collaboration 
obtain a value for the inverse lattice spacing, a -1 
= 2.4(1) GeV. The Q value for the fit is good at 
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Figure 1. The spectrum for bb states plotted 
relative to the mass of the T (IS). Solid hori- 
zontal lines mark the experimental data for the 
3 Si and for the spin-average of the 3 Po,i,2 (Xfe) 
ground states and radial excitations. Grey hori- 
zontal lines mark the expected masses of the 1 P± 
(hb) states. The vertical scale is in GeV. Squares 
mark the results from the NRQCD collaboration. 
The stars are those from the UKQCD collabora- 
tion (using leading order NRQCD). Their results 
are all for spin-averaged masses. The cross on 
T(1S) shows that its mass was assumed in the 
results. 



0.4. Taking the ratio of plaquette values at (3 = 
6.0 and (3 = 6.2 and using lattice perturbation 
theory gives a value for the inverse lattice spac- 
ing at = 6.2 of 3.2(2) GeV. It is this value of 
a -1 that has been used to convert the UKQCD 
results from lattice to physical units. The FNAL 
group have calculated so far a value only for the 
IP- IS splitting for bb and they use this to deter- 
mine a -1 , giving no independent predictions to 
appear in Figure 1. They obtain 2.7(2)GeV at f3 
= 6.1, a value in agreement (using perturbative 
extrapolations) with the NRQCD value above. 

Notice that all these values of a -1 are larger 
than those quoted at similar /3 values from the 
string tension or light hadron spectra. This sim- 
ply reflects the scale dependence of the determi- 



nation of a -1 from the quenched approximation, 
when the long distance heavy quark potential is 
too steep and the short distance potential not 
steep enough. The difference is that the short 
distance potential can be corrected perturbatively 
for its quenching errors. 

Figure 1 shows a good fit to experimental 
data [jl3| within the errors. The statistical and 
systematic errors of the NRQCD calculation || 
are such that it is almost possible to see the effect 
of the quenched approximation. The IP level is 
slightly too low and the 2S slightly too high. Es- 
timates from a lattice Schrodinger equation de- 
scribed above give a correction for quenching and 
gluonic 0(a 2 ) effects that produces a very good 
fit to experiment. 

Figure 2 shows the hyperfine spectrum from the 
NRQCD collaboration for bb IP states compared 
to experiment. Separate results are given for dif- 
ferent components of the 3 p2 and 3 Pi states. The 
fit is good, and provides a stringent test of the 
coefficent of the a.(D x E) term in the spin- 
dependent action. Using tree-level coefficients 
without correction for tadpoles would give split- 
tings reduced by a factor of 50%. Values for 
the (unobserved) IS hyperfine splitting are 29(4) 
MeV for NRQCD, 25(4) MeV for FNAL and 18(1) 
MeV for UKQCD using heavy Wilson with a tree- 
lvel improvement coefficient. It seems likely that 
the UKQCD value is too low. 

A pole mass can be extracted for the b quark 
from a perturbative analysis of the dispersion re- 
lation for the T at finite momentum. A fit to 
the energy of propagators at zero and non-zero 
p gives the T mass in lattice units and its zero- 
momentum energy (Mi and Mi in equation 7). 
However, M\ and M2 are should be related by 

M 2 = 2(Z M M° -E ) + Mi (7) 

where Zm and E are perturbatively calculable 
mass renormalisation and energy shift parame- 
ters and is the bare heavy quark mass. The 
NRQCD collaboration § compares the results for 
M2 from equation 8 and the value extracted di- 
rectly using equation 7. They agree to 1% for a 
variety of bare heavy quark masses. This agree- 
ment is at the level expected for a nlo calculation 
and should be seen as a triumph for the whole 
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Figure 2. The hyperfine spectrum of the bb IP 
states, relative to the spin-average of the 3 -Po,i,2 
(Xb) states. The horizontal lines mark experimen- 
tal data; the vertical scale is in MeV. The squares 
are results from the NRQCD collaboration, with 
error bars where these would be visible. 



method. It requires both that the nlo coefficients 
be correct and that the lattice perturbation the- 
ory is well-behaved. Thus the assumption of tad- 
pole dominance of the coefficients is correct and 
there are no unexpectedly large perturbative or 
non-perturbative contributions to Zm or Eq. The 
numbers used for Zm and Eq have been calcu- 
lated to 0{g 2 ) for the action used by the NRQCD 
collaboration |1(J and using the scheme and scale 
for g 2 advocated by Mackenzie and Lepage ||. 

Fixing M° so that the experimental value for 
Mr is obtained then yields a value for the renor- 
maliscd or pole b mass, ZmM®. This value is 
4.7(1) GeV where the main source of error is that 
of higher order corrections to Eq . Simulation er- 
rors are not significant. 

2.2. cc results 

New results for the cc spectrum were presented 
only by the NRQCD collaboration. Their results 
are shown in Figure 1 of the talk by A.J. Lid- 
sey [Q. Now experimental results are available 
for both the rjc^So) and the h c ( Pi). The spin- 



averaged spectum is therefore presented with the 
spin-averaged s state at zero. There are two s 
states below threshold for decay to heavy-light 
channels. The 2S state is fit using a 3-exponential 
fit to 2 correlation functions. The fit to experi- 
ment is good and certainly within systematic er- 
rors of 10%. An estimate for the position of a 
d state was obtained by calculating a correlation 
function for the 1 Z?2- This is compared to the 
posited 3 L>i state, -0(3770). The 1 D 2 is higher, 
as expected. There may also be significant cou- 
pling to decay channels here which would tend to 
distort the spectrum in the quenched approxima- 
tion. 

Figure 3 shows the cc hyperfine spectrum. 
Both s and p states are given relative to their 
respective spin-averages. For the s state hy- 
perfine splitting there is some disagreement be- 
tween groups. The NRQCD collaboration ob- 
tains a value of 96(4)MeV, the FNAL collabora- 
tion, 64(4)MeV and the UKQCD collaboration, 
50(l)MeV. The experimental value is 116MeV. 
There is a 30% systematic error in the calcula- 
tions from higher order terms in the quark action 
and also a systematic effect from quenching which 
tends to reduce the value. If the effect of quench- 
ing is a 40% reduction then the NRQCD value 
looks high. If a 20% effect then the other val- 
ues look low. The NRQCD result has the small- 
est systematic error from the quark action, which 
should give the best accuracy in determining the 
quark mass (something the hyperfine splitting is 
sensitive to). On the other hand the FNAL col- 
laboration have results for the hyperfine splitting 
for different values of (3 and see no strong system- 
atic effect from 0(a 2 ) terms that they have not 
included. This disagreement needs to be resolved. 

3. CONCLUSIONS 

Heavy quark spectroscopy can provide a strin- 
gent test of QCD. The bb spectrum can be calcu- 
lated to the 1% level using an action correct to 
nlo in v 2 /c 2 . Such an action can be provided 
within an NRQCD or heavy Wilson approach. 
Current results in NRQCD are at this level in the 
quark action. There are remaining systematic er- 
rors present from the gluon configurations used. 
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Figure 3. The hyperfme spectrum of the cc IS 
and IP (x c ) states, relative to their respective 
spin-averages. The horizontal lines mark exper- 
imental data; the vertical scale is in MeV. The 
squares are results from the NRQCD collabora- 
tion, the circles from the FNAL collaboration and 
stars from the UKQCD collaboration (heavy Wil- 
son). Error bars are shown where they would be 
visible. 



These should be corrected in future calculations. 
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